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ABSTRACT 

Aims. We present the results obtained from a long-term spectroscopic campaign devoted to the multiplicity of O-type stars in the 
young open cluster NGC 2244 and in the Mon 0B2 association. 

Methods. Our spectroscopic monitoring was performed over several years, allowing us to probe different time-scales. For each star, 
several spectral diagnostic tools are applied, in order to search for line shifts and profile variations. We also measure the projected 
rotational velocity and revisit the spectral classification. 

Results. In our sample, several stars were previously considered as spectroscopic binaries, though only a few scattered observations 
were available. Our results now reveal a more complex situation. Our study identifies two new spectroscopic binaries (HD46149 in 
NGC 2244 and HD 46573 in Mon OB2). The first object is a long-period double-lined spectroscopic binary, though the exact value of 
its period remains uncertain and the second object is classified as an SB 1 system with a period of about 10.67 days but the time series 
of our observations do not enable us to derive a unique orbital solution for this system. We also classify another star as variable in 
radial velocity (HD 46150) and we detect line profile variations in two rapid rotators (HD 46056 and HD 46485). 
Conclusions. This spectroscopic investigation places a firm lower limit (17%) on the binary fraction of O-stars in NGC 2244 and 
reveals the lack of short-period 0+OB systems in this cluster In addition, a comparison of these new results with two other well- 
studied clusters (NGC 6231 and IC 1805) puts forward possible hints of a relation between stellar density and binarity, which could 
provide constraints on the theories about the formation and early evolution of hot stars. 

Key words, binaries: spectroscopic - stars: early-type - open clusters and association: individual: NGC 2244 - open clusters and 
association: individual: Mon OB2 



1. Introduction 

O-type stars are preferentially found in young open clusters or 
in OB associations, suggesting these two latter environments to 
play a key role in massive star formation. Detailed observations 
could provide interesting constraints enabling to choose between 
the different formation scen arii proposed by theoret icians (for 
a summary of these, see e.g. IZinnecker & York3l2007l) . To this 
aim, massive stars in clusters are a target of choice because they 
represent a homogeneous population (same distance, same age 
and same chemical composition). Moreover, the study of the 
population of O-type stars in open clusters, and in particular of 
their multiplicity, is of major interest for the understanding of 
the physical properties of these objects. Other parameters (or- 
bital period, eccentricity, mass ratio,...) of such multiple systems 
but also the binary frequency can provide a unique signature of 
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star formation and dynamical interactions within the cluster dur- 
ing the very early phases of its existence. 

Investigations based upon inten se spectroscopic c ampaigns 
in some clusters (e.g . NGC 6231, SanaetaD I^QS'; IC1805, 
iDe Becker et al.ll2006l) have been performed over several years. 
These detailed obs ervations lead to a serious c orrection of the 
results proposed bv lGarcfa & Mermilliodl (12001}) for O-type star 
rich clusters, by greatly improving our knowledge of the precise 
physical parameters of the binaries in these clusters. In order to 
extend these studies, we have undertaken a similar spectroscopic 
monitoring of the O-type star population in another young open 
cluster, NGC 2244. 

With an inferred age of about 2-3 Myr jChen et all 

I2007h and a radial v elocity (RV) of 26.2 + 3.4 kms"' 



dKharchenko et all |2005|) . NGC 2244 lies inside the Rosetta 



nebula and forms the core of the Monoceros OB2 association. 
This cluster is the youngest of two or thre e subgroups of OB stars 
and stellar aggregates in this association dLi & Smithll2005l) . 

Some stars in NGC 2244 were observe d by 
iPlaskett & Pearcd (|1931|) and reanalyzed by Petrie & Pearcd 



(1196 Ih . These authors found that the majority of these ob- 



jects has quite constant RVs. These results were, however. 
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subseq uently challenged by lAbt & Biggsl ( Il972h and iLiu et alj 
([19891). 

From a compilation of heterogeneous data, 
iGarcia & MermiUiodI (1200 1[) estimated the proportion of 
spectroscopic binaries in NGC 2244 to be about 50% of the 
O-type stars (three out of six considered as binaries). Such 
an O-star binary fraction would make this cluster the third 
binary-rich cluster in their sample after NGC 6231 and IC 1805. 
As already mentioned above, the fi gures for these two clusters 
were significant l y revised recently ( Rauvv & De BeckeA 120041: 
iDe Becker etalL l2006t ISana et al.L l2008h . In addition, up to 
now, only two spectroscopic binaries have been d etected with 
certainty in Mon OB2 : HP 47129 (P laskett's star, iLinder et all 
120081) and HD 48099 (ISticklandlll996 ; Mahy et al. in prepara- 
tion) and both of them lie outside NGC 2244. These conflicting 
results clearly call for a new, thorough analysis of the binarity in 
NGC 2244. 

Other studies in different wavelength domains reveal addi- 
tional information on the content and overall organization of 
the cluster. Chandra mosaic observations of NGC 2244 were 
presented by Wang et al. (2008). The X-ray spectra of the OB- 
type stars were found to be rather soft and consistent with 
the standard model of small-scale shocks produced, in the in- 
ner regions of the s t ellar winds, by the line-driving instabil- 
ity dFeldmeier et al.L Il997h . The X-ray luminosities of these 
stars follow the canonical Lx/^^boi relation (Sana et al., ,2006b) 
rather closely and reveal no evidence for extra emission due 
to wind-wind collisions in binary systems. I Wang et akl (l2008l) 
found that the fainter X-ray sources are strongly concentrated 
around HD 46150 (05.5 V) whilst no strong clustering of X- 
ray sources was found around the most massive component 
HD 46223 (04 V). Following these authors, there are two possi- 
ble explanations for this difference: either HD 46223 was ejected 
as a result of dynamical interactions (but it does not display a 
strong proper motion) or it may actually be younger and would 
not be part of the same population as the central part of the clus- 
ter. 

The present paper, based on a long-term spectroscopic cam- 
paign, is organized as follows. Section 2 describes the observ- 
ing campaign. The properties of the individual objects (spectral 
classification, projected rotational velocity and multiplicity) are 
presented in Section 3. The multiplicity and the rotational ve- 
locities of massive stars in NGC 2244 are discussed in Section 4 
while Section 5 provides a summary of our results. 

2. Observations and data reduction 

Our team collected 207 spectra of 9 stars over 9 years using sev- 
eral different telescopes. This sample allows us to investigate the 
presence of short and long-term variations for each star. 

Most of our observations were obtained with the 1 .52m tele- 
scope of the Observatoire de Haute-Provence (OH?) equipped 
with the Aurelie spectrograph. The detector used was a thin 
back-illuminated CCD with 2048 x 1024 pixels of 13.5 yum^. 
The grating #3 (600 lines/mm) allowed us to obtain spectra with 
a resolving power of 10000, centered on 5700 A and cover- 
ing wavelengths between 5500 A and 5920 A. Exposure times 
were typically 10 to 45 minutes resulting in signal-to-noise ra- 
tios larger than 150. Such high quality observations were needed 
to search for weak variations or for the spectroscopic signature 
of faint companions. We also collected spectra with the same 
instrument in the blue domain (centered at 4670 A and cov- 
ering the 4450-4900 A region). These data were spread over 



28 days in November 2007 and over 6 days in September and 
October 2008. The whole dataset was reduced in a classical way 
(bias and flatfield corrections) using the MIDAS software de- 
veloped at ESO (European Southern Observatory). The spectra 
were wavelength calibrated using Thorium-Argon comparison 
spectra obtained immediately before or after the stellar spec- 
trum. The spectra were finally normalized to the continuum by 
fitting polynomials of degree 4 or 5. 

Additional spectra were obtained at Observatorio 
Astronomico Nacional of San Pedro Martir (SPM), in Mexico, 
with the 2.1m telescope equipped with Espresso. This echelle 
spectrograph gives 27 orders in the 3780-6950 A wavelength 
domain with a resolving power of R = 18000. The CCD detector 
was a SITE 3 optical chip with 1024 x 1024 pixels of 24 fim^-. 
Typical exposure times range from 5 to 15 minutes. We had to 
add consecutive spectra of a given night at the expense of time 
resolution to obtain signal-to-noise ratios close to 300. The data 
were reduced using the echelle package available within the 
MIDAS software. 

In addition, we gathered spectra from both Elodie and ESO 
archives. Elodie was an echelle spectrograph giving 67 orders 
in the 3850-6850 A spectral range (R = 42000) mounted at the 
1.93m telescope at OHP. We retrieved 13 spectra spread be- 
tween November 1999 and November 2005. The mean signal-to- 
noise ratio for our targets was 150 for an exposure time ranging 
from 10 min to 1 hour. The ESO archives provided us with 11 
spectra taken with FEROS (Fiber-fed Extended Range Optical 
Spectrograph) mounted at the ESO/MPG 2.20m telescope at La 
Silla (R = 48000). The data reduction was performed using a 
modified FEROS pipeline working under the MID AS environ- 
ment. B eyond the modifications already described in lSana et al.l 
(I2006ah . several new features were implemented. We used a 2D 
fit of the order position for improved stability. We also used the 
wavelength calibration frames obtained with the "new" ThArNe 
calibration lamp (see e.g. FEROS-II user manual 77.0, §2.5). 
The latter however heavily saturates in the red, even on the short- 
est exposures planned by the calibration plarQ. Therefore, we 
computed a master calibration frame for each exposure time. 
We then proceeded with the extraction of the calibration spectra 
separately for each master frame. We finally used the detected 
lines from the 3s/15s/30s master wavelength calibration frames 
to respectively calibrate the orders 30-39/10-29/1-9. The satu- 
rated lines or the lines outside the linearity range of the detector 
were rejected from the fit and the detection threshold was op- 
timised for each order. Finally, a 2-iteration 3-sigma clipping 
method was used to discard the very few remaining discordant 
lines (mostly because of artifacts related to poorly corrected bad 
column effects). The obtained wavelength calibration residuals 
were all in the range 2.7-2.9 mA. 

At last, we retrieved two spectra from Asiago archives, 
with a mean signal-to-noise ratio of about 155, to increase the 
dataset of HD 46149. These two spectra were taken by the 1.82m 
Copernicus telescope equipped with the echelle instrument (R = 
28600). Unfortunately, no flatfields were available and we could 
only apply an overscan correction. Due to the lack of a well- 
suited calibration atlas, the spectra were approximatively cali- 
brated in wavelength but this first order calibration was subse- 



' As a reminder, the latter provides 3 series of two exposures with 
increasing exposure times of 3s, 15s and 30s, a scheme that is repeated 
at the beginning and at the end of the night. 
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quently checked using the interstellar lines. The overview of the 
observations is given in Table [lIFI (available electronically). 

For most stars in NGC 2244, the radial velocities (RVs) 
and the equivalent widths (EWs) were determined by fitting 
Gaussians on the line profiles, except for the two rapid rotators. 
Indeed, Gaussians properly fitted the observed profile of moder- 
ately rotationally broadened lines in massive stars. To measure 
the RVs, we fitted the bottom half of the line profiles while the 
entire line was used to estimate the EWs. For the two rapid rota- 
tors, we used a synthetic line profile generated with the projected 
rotational velocity of the star to determine the Doppler shifts and 
we integrated the complete line profiles to measur e the EWs. 
We no te that we used the same rest wavelengths as IConti et aTl 
(1 19771) to calculate the RVs. 

Moreover, we inferred the spectral type and the luminos- 
ity class of each star by using the classification criteria of 



ity class 01 each s tar by usmg t he cla ssiiica tion crit e ria or 
Conti & AlschuieJ [l9m . IContil (IT973h and iMathvsl ([T9881 



19891) hereafter referred to as "Conti's criteria". To support 



our results, we visual l y com pared the spectra to the atlas of 
IWalborn & Fitzpatrickl (1199 0'). We also checked the luminos- 
ity class by comparing the computed visual absolute magnit ude 
with the photometric calibration proposed by IMartins & Pled 
(12006.). To this aim, we need a value for the distance of the clus- 



ter. 



Massev et al.l d 19951) derived a spectroscopic parallax to this 
region of 1.9+0.1 kpc, more recent photometric studies place 
NGC 2244 at a distance of 1 .4 to 1 .7 kpc ( Hen sberge et al., 2000) 



in agreemen t with previous estimates (lOgura & Ishidal 119811: 
lPerezl[T99Th . In the present paper, we adopt a mean distance for 
NGC 2244 of 1.55 kpc. However, we note that a change of dis- 
tance by 10% do not modify our conclusions. Finally, we com- 
pared our spectral classification with previous re sults from e ar- 
lier studies (notably quoted in the Reed catalo gue lReedll2005l) . 



3. The 0-type star population in NGC 2244 and l\1on 
0B2 

Our study of the multipUcity of O-type stars is based on two 
wavelength domains: a blue region (4450^900 A) and a yel- 
low one (5500-5920 A). In addition to the RV and profile anal- 
ysis, the former setting provides an accurate determination of 
the spectral classification of each star; the numerous interstellar 
features of the latter setting enable to check the quality of the 
wavelength calibration and RV measurements. We provide, in 
Table 2 (available electronically at CDS), the heliocentric Julian 
days (HJD) and the RVs of each spectrum while Table[3]quotes, 
object by object, the mean RVs and the 1-cr standard deviation 
of the different lines. The data are expressed inkms ^ 

In the following, an object will be considered as a true spec- 
troscopic binary (Section 13. Il l if we detect either the presence 
of significant, periodic RV variations or the signature of a com- 
panion whose spectral lines move in anti-phase with those of the 
main star. If a star fulfills the former criterion, a full orbital solu- 
tion can be computed, leading to the precise determination of the 
orbital parameters. When the latter criterion is fulfilled, it only 
allows us to evaluate the spectral type of the two components and 
their mass ratio, with only loose constraints on the orbital period. 
If none of these criteria is fulfilled, but significant RV variations 



^ The first and the second columns give the name of the object and the 
membership to the cluster or association (Ogura & Ishida, 1981). The 
next columns are the epochs of the campaign, the telescope/instrument 
used, the observed wavelength domain, the number of spectra obtained 
and the time elapsed between the first and the last spectrum of the run, 
expressed in days. 



are detected, then the star will be considered as a "RV variable", 
and thus a binary candidate. The RV variations are considered 
as significant if they are larger or equal to three times the av- 
erage RV measurement error. For all stars except the two rapid 
rotators, this error is well represented (as was confirmed by tests 
on individual data) by the RV dispersion of the narrow diffuse 
interstellar band (DIB) at A 5780 which has similar strength and 
width than the stellar lines. For the rapid rotators, evaluation 
from shifts of the synthetic line profile showed the RV error to 
be about 10 kms"'. Object for which significant RV variations 
are detected is presented in Section [T2l The remaining objects 
are presented in Section 1331 

3. 1 . Spectroscopic Binaries 
3.1.1. HD46149 

The first spectro s copic study of HD46149 was done by 
Pla skett & Reared (1193 lb who measured a mean RV of about 
45.0+ 1 .5 km s~ ' . iPetrie & Reared (1 1 96 ll) confirmed, on the ba- 
sis of new data, that the star might be single. They obtained a 
mean RV of about 3 8 . ± 2 . 1 km s~ ' . iLiu e t al] d 1 9 89 ) analyzed 
three datasets obtained over three different epochs (November 
1987, February 1988 and October 1988). They found a RV dis- 
persion of about 10.2 km s ' and, t herefore, considere d the star 
as a potential binary. More recentlv. fTurner et"aLl (|2008|) reported 
some velocity variabilities of HD 46149 due to the presence of a 
close spectroscopic companion. 

In February 2006, the measurements of RVs for six spectra 
over six consecutive nights did not show any significant variation 
(RV =21.6 + 2.1kms '). Moreover, no variation was detected 
in the line profiles. However, if we compare these RVs with the 
April 2007 data (HJD ~ 2 454 199, in Fig.[I]i, we obtain a shift 
close to 25 km s ' for the He i A 5876 line, and a RV shift of about 
40 km s"' for He i A 447 1 is measured between the data taken in 
January 2006 and in November 2007. Our analysis of the whole 
spectroscopic campaign on HD 46149 clearly shows significant 
RV shifts between the observing runs. The two extreme RV val- 
ues in our sampling are approximately -5 (Jan. 2006, Mar. 2008) 
and +50kms-i (Apr. 2007, Nov. 2007 and Sept.-Oct. 2008). In 
addition, the FEROS spectrum and the data collected between 
January and April 2008 allowed us to detect the presence of 
a companion in the red wing of most of the He: lines (visi- 
ble in Fig.[I]at, for example, HJD ~ 2453 739 and HJD ~ 
2454545). The determination of the individual RVs enables 
us, in this case, to estimate the mass ratio of the two compo- 
nents. At the largest observed separation (at HJD ~ 2454545), 
the best precision on the RVs of both stars was obtained from 
analyzing the He 1/15876 line, which yields velocities of -5.8 
and +143.3 km s"' for the primary and the secondary, respec- 
tively. The RV of the blended spectrum, +52.5 km s ' at HJD ~ 
2452 664, was measured on the He 1/15876 line. From these 
shifts, we estimate a mass ratio (M1/M2) of about 1.6. We nev- 
ertheless emphasize that the phases for which we measured the 
RVs might not be the extreme ones of the orbital cycle. In con- 
sequence, our mass ratio computation must be considered pre- 
liminary at best. We also estimated the flux ratio between the 
primary and the secondary. We measured, from the deblended 
He I /1 5876 line, an EW for the primary of 572 mA which we 
compared with the average EW quoted by Conti (1973) for such 
stars. We obtained an optical flux ratio L\I(L\ + L2) for the pri- 
mary of about 0.7. In consequence, the flux ratio L1/L2 is close 
to 2.3. The obtained mass and flux ratios thus suggest an early 
B-type for the secondary (for comparison, from the stellar pa- 
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Table 3. The mean radial velocities and the 1-cr dispersions. 



Star Hei^4471 HeiiA4542 HeiU4686 Oiii.t5592 DIBi5780 Civi5801 Civ^5812 Hei/i5876 Naii5890 

HD46056 36.1 ±6.9 49.6 + 9.3 30.4±12.1 ... 26.0±4.1 ... ... 28.7±9.7 23.3±2.0 

HD46149 27.3±21.8 27.6± 18.7 28.4±22.0 12.6±214 23.5±4.7 24.1 ±21.3 20.7±20.3 22.3±194 214±3.0 

HD46150 31.2±4.0 45.3±8.1 44.8±5.0 26.9±7.8 23.8±3.9 48. 8± 11.7 34.0± 12.9 32.6±54 20.7±24 

HD46202 39.8±3.1 38.8±3.8 38.0±3.2 31.8± 1.3 25.4±1.5 404±2.2 36.1 ±2.2 40.0±14 23.5±1.2 

HD46223 34.5±6.7 44.5±3.2 47.0±3.3 31.1±74 23.8±4.6 50.9±2.9 38.6±64 35.6±2.7 234±3.4 

HD46485 31.5±7.5 43.5± 10.7 34.0±94 ... 21.8±5.5 ... ... 30.8±11.1 21.3±4.0 

HD46573 47.0±6.6 48.4± 10.0 50.2±7.6 47.9±8.5 27.6±3.9 59.5±5.7 53.2±6.1 50.3±54 26.9±1.4 

HD 46966 39.2 ±14 44.9 ±2.0 45.3 ±2.0 384 ±2.6 23.4 ±4.7 47.7 ± 3.0 434 ±3.3 40.8 ±24 20.8 ±2.2 

HD48279 314±3.6 40.6±2.6 38.9±5.1 32.8± 11.8 23.1 ±5.0 43.3±6.0 31.9±74 39.5±2.1 20.8±3.3 



rameters listed in lMartins et al.ll2005l a binary system composed 
of an 08 V and an 09.5 V would have a theoretical mass ratio of 
M1/M2 =1.3 and a flux ratio L1/L2 of about 1.9). 




4470 4480 5870 5880 

A (A) X (A) 



Fig.l. He I i 4471 and He 1/15876 line profiles of HD 46149 ob- 
served during our spectroscopic campaigns. The time is given 
on the right side of the figure in the format HJD - 2 450 000. 
All spectra are displayed in the heliocentric frame of reference. 
The vertical dashed line represents the rest wavelength. The sig- 
nature of the secondary is clearly seen, in particular at HJD ~ 2 
454 545. 



varies according to the studied line. For example, the computed 
periodogram corresponding to He 1/15876 yields a period of 
about 47 days while the one corresponding to Civ /1 58 12 in- 
dicates a 1600 day period. We also performed a period search 
on the combined RV list of He i /1 447 1 and He i A 5876. Even in 
this case, the detected peak in the periodogram is not outstanding 
enough to support a particular orbital period. The data sampling 
is thus too poor at the present time to constrain the period of 
HD 46149. 

In order to separate the two components of HD 46149, we 
used our disentangUng prog ram, based on the method pro- 
posed bv iGonzalez & Levatd (12006). This consists in an itera- 
tive process which alternately uses the spectrum of one compo- 
nent (shifted according to its radial velocity) to remove it from 
the observed spectra and to calculate a mean spectrum of the 
other component. Furthermore, this technique allows us to com- 
pute the RVs of each star by cross-correlation even at phases 
for which the lines are heavily blended. We applied the dis- 
entangling algorithm in the wavelength domains 4000^220 A 
and 4450^730 A. To create the RV cross-correlation mask, 
we used the He U 4026, Si iv /1 4089, He 1/14120, He 1^4143, 
He I A 447 1 , Mg n /1 448 1 , He i /1 47 1 3 Unes for the secondary and 
we added the Si iv /1 4 1 1 6, He 11 /1 4200, He 11 A 4542 fines for flie 
primary. The calculated mean spectrum o f the secondary star 
looks similar to that of a BO V star in the Walb orn & Fitzpatrickl 
(Il990i) atlas. However, because of a poor sampling of the data, 
this result must be considered very preliminary. 

HD46149 was classified as an 08.5V((f)) star 
(iMassev et al.L 1 19951) but we revise the classification to 08 
V. Note that this does not depend on the degree of blending of 
the spectra. With the mass ratio computed above and a theoret- 
ical mass for the primary star close to 22 Mq dMartins et al.L 
i2005i) . we estimate the secondary mass at approximately 14 
Mq. Therefore, the system probably consists of 08 V -1- BO-1 
V stars: HD46149 is the first O+B binary discovered so far in 
NGC 2244. 



We then applied the Fourier analysis on the time se- 
ries of primary RVs of He 1/14471, He 11 /1 4542, CiV/15812 
and H e I /1 5876 followi ng the meth od described bv iHeck et al.l 
( Il985h and revised bv iGosset et al.1 fcOOl). In all cases, we ob- 
served the highest peaks in the periodograms but its position 



3.1.2. HD 46573 

This star i s a member o f Mon OB2, but does not belong to 
NGC 2244 dTurneii [19761) . According to the RV variability crite- 
rion, the RV dispersions for the stellar lines of HD 46573 are at 
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the limit of being considered as significant. Ho wever, we com- 
puted the Temporal Variance Spectrum (TVS, iFuUerton et alj 
1 19961) from the Aurelie spectra for the He i A 447 1, Hen A 4542, 
C IV /l/l 5801-12 and He 1/15876 lines. We obtained, in each 
case, double peaked profiles (Fig|2]i. Together, the RV dispersion 
and the double peaked TVS structures lend signifi cant support 
to the binary nature of the star, first envisaged by iMason et alj 
(119981) . 
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Fig. 2. The mean spectrum and TVS of HD 46573 com- 
puted from the Aurelie spectra for the He 1/14471 (top) and 
He 11/14542 lines (bottom). In both cases, the double peaked 
structure is clearly visible. The dotted line illustrates the 99% 
significan ce level for the variab ility evaluated following the ap- 
proach of iFuUerton et al.1 (Il996l) . 



We performed a Fourier analysis of the RVs. We indepen- 
dently applied this method to three RV time series (He 1/14471, 
He 11 A 4542 and He i A 5876) but the semi-amplitude spectra did 
not allow us to identify any dominant peak likely to be asso- 



ciated with a well-determined period. We then combined the 
RVs of the He 1/14471 and He 1/15876 lines. The highest peak 
in the resulting semi-amplitude spectrum occurs at a frequency 
of 0.0937 d"' (Fig.O which corresponds to a period of about 
10.67 days but this value is still very preliminary since ran- 
dom fluctuations can not be totally excluded. The SBl orbital 
solution of HD 46573 is calculated with the program LOSP 
(Liege Orbital Solution Package, Sana&Goss et, A&A submit- 
ted) based on the method of Wolf e et alj (Il967l) . We assigned the 
same weight, i.e. 1.0, to all spectra and fixed the period value to 
10.67 days. Table|4]gives the main orbital elements computed for 
HD 46573. To (expressed in HJD-2450000) refers to the time 
of periastron passage for e + ^ and to the conjunction with the 
primary star in front for e-^.y,K and a sin / are respectively 
the systemic velocity, the semi-amplitude of the radial velocity 
curve, and the projected separation between the centre of the star 
and the centre of mass of the binary system. The orbital elements 
are given respectively for a circular and for an eccentric orbit 
from flie RVs of bofli He 1/14471 and /1 5876. Folded with this 
period, our data favored a non-zero eccentricity (Fig.|4l bottom 
panel). While not in sufficient number to calculate a full orbital 
solution, the RVs from Hen and Civ lines support the results 
from the Hei lines (Fig.|4|. However, there are still some phase 
intervals without observations and we therefore need more data 
to secure the period and the orbital parameters. 
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Fig. 3. The square root of the power spectrum (top) and the spec- 
tral window (bottom) of HD 46573 computed from RVs of both 
He I /1 4471 and /1 5876 lines. The y-axis corresponds to the semi- 
amplitude of the signal. The frequency is expressed in d ' . 



We then computed , with the Fourier transform method 
jSimon-Dfaz & Herrerol |2007|) . the projected rotational veloc- 
ity of HD 4657 3. Prev iously estimated at about 95kms"' 
(IConti & Ebbeti. [T977h . we measured, for four lines 
(Heii/14686, OiiU5592, CivA5801 and Hei/15876), a 
mean value of about v sin ; =110+ 18km s"' . 
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Fig. 4. Top: the best-fit circular radial velocity curve for 
HD 46573 computed from the RVs of the He 1/14471 (open 
squares) and /1 5876 (filled squares) lines for an orbital period 
of 10.67 days. The open triangles correspond to the RVs of 
the He 11/14542 line and the stars to the RVs of the Civ /1 58 12 
line. The eiTor bar corresponds to ± 1-cr RV dispersion of the 
DIB A 5780. Bottom: the best-fit eccentric orbital solution for 
HD 46573 computed for the same two lines. 

Table 4. Possible orbital parameters of the SBl solution of 
HD 46573. 



Parameter 


Circular 


Eccentric 


P(days) 


10.67 (fixed) 


10.67 (fixed) 


e 




0.47+0.13 


ij (degrees) 




254.63+14.2 


To 


3732.24±0.26 


3731.68+0.48 


y (km s"') 


49.7+ 0.7 


50.9+0.8 


^:(km s-') 


6.9±0.9 


8.5+1.1 


a sin i (Rq) 


1.46+0.19 


1.58+ 0.24 


f(m) (Mo) 


0.0003+ 0.0001 


0.0005+0.0002 


r.m.s. (km s"') 


3.5 


2.6 


ProbO(^2 >) 


78% 


99% 



The spectr al type was previ ously quoted as 07.5 V ((f)) 



(IConti&Leed Il974j). 
OlllKif)) (IWalbornL " 



07V (iBisiacchi et al.L IT982h 



19711) . In our data, the Hei/i4471 to 



He 11 /1 4542 ratio impUes an 07.5 type while the EW ratio 
between Si iv /1 4088 and He 1/14143 favors a Main-Sequence 
(MS) luminosity class. The visual absolute magnitude is 
more uncertain since the distance to the star is not exactly 
known. Assuming that Mon OB2 lies at the same distance 
as NGC 2244, we estimated My for the two recent extreme 
values of the distance of NGC 2244 (1.4 kpc and 1.7 kpc). We 
obtain My = -4.6 7 and My = -5.07, using V = 7.933 and 
E{B - y) = 0.61 (iMafz-Apellaniz et al.L |2004|). The reported 
result s are closer to a MS star than to a giant (Martins & Plezl 
|2006|) . The spectrum also shows weak Niii/1/1 4634-41 emis- 
sion lines and a strong He n A 4686 absorption. In Ha, we detect 
a weak emission in the absorption line probably due to the 



surrounding nebula. Thus, we favor an 07.5V((f)) spectral 
type, for HD 46573. 

In summary, it exists two possible explanations for which the 
secondary component is not detected: a low-mass companion or 
a small inclination of the orbital plan. However, with such a large 
V sin / for the primary, it seems to be unlikely that the system is 
seen under a low inclination angle. In consequence, we assume 
that HD 46573 is an 07.5 V((f)) star which possesses a lower 
mass companion. 

3.2. Radial velocity variables 
3.2.1. HD46150 

Reported as the second hottest star in NGC 2244, HD46150 
is also one of the most st udied objects in t h is clu ster. This 
star was first observed by iPlaskett & Reared (1193 lb who re- 
ported constant RVs. Since the 1 970's, many other inv estigations 
were devoted to HD 4 6150 (Abt & Biggsl [T972HConti & Ebbetsl 
119771; iGarmanv et al.L [l980) but the interpretations of the mul- 
tiplicity of the star g reatly varied from one author to another. 
iGarmanv et al.l (Il980l) considered the changes in the RVs to be 
a cons eque nce of motions in the stel lar atmosphere. iLiu et al.l 
( Il989h and lUnderhill & Gilrovl (11990) suggested that the star 
could be a spectroscopic binary even though the sign ature of 
the se condary component had not been detected yet. F ullertonI 
(Il990h observed asymmetric line profiles, typical of an SB2. 

Between November 1999 and October 2008, we have col- 
lected 34 spectra spread over different time-scales. We ob- 
served a slight RV variability in the absorptions (Fig.|5] no- 
tably compares the Hen /1 4542 line at HJD~2453 772 and 
HJD ~ 2 454 397) as well as in the emission lines which is 
supported by the larger dispersion for stellar lines compared 
to the interstellar ones. For instance, we noted a RV shift of 
He II A 4542 in the blue data by 28 km s"' and of He i /1 5876 in 
the yellow spectra by 21 km s Moreover, the RV dispersion 
is smaller for the He i lines than for the He ii and metallic lines 
which are clearly RV variable. The presence of an unresolved 
companion could explain this effect. Indeed, if the secondary is 
cooler, the Hei lines of the primary will be more affected by 
blends with the lines of the secondary component. The RV varia- 
tions of lines associated with a stronger ionization (He ii, C iv,...) 
will thus reflect better the true motion of the primary which leads 
to a larger RV dispersion. 

In order to compute the TVS from all spectra of HD 46150, 
we had to degrade the spectral resolution of the echelle spec- 
tra down to that of the Aurelie data. We proceeded by convolv- 
ing the echelle spectra with Gaussians to obtain a homogeneous 
spectral resolution for all spectra. Variations in double peaks 
(in the Hei/14471, Heii/14542, He ii /1 4686 and Civ /l/l 5801- 
12 lines) and even variations in triple peaks, shown in Fig.|6] for 
the He I /1 5876 line, are observed. Such features are reminiscent 
of what is usually found in spectroscopic binaries. However, no 
obvious signature of a companion could be found. Moreover, we 
note that these variations in double peaks are not visible on short 
time-scales (~ 7 days), suggesting a rather long-period binary. 

To check for the presence of a putative period, we applied the 
Fourier method to the measured RVs. We independently used 
this technique for the Hei/14471, Hen/14542, CiV/15812 and 
He I /1 5876 lines. For each line, several peaks are present in the 
periodogram but their semi-amplitude nearly corresponds to the 
standard deviation for the RVs of the interstellar lines. We also 
computed the power spectrum by combining the RVs of the 



L. Mahy et al.: Early-type stars in the young open cluster NGC 2244 and in the Mon 0B2 association 



7 





4540 4550 
A (A) 



5880 
K (A) 



Fig, 5, He n ^ 4542 and He i /( 5876 line profiles of HD 46150 ob- 
served during our spectroscopic campaigns. 



He I /i 4471 and the He 1/15876 lines but, once again, the peak 
is not significant enough to ascertain the period. 
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Fig. 6. Mean spectrum and TVS of HD46150 for the 
C IV AA 5801-12 and He i A 5876 Unes. The dotted line illustrates 
the 99% signifi cance level for the vari ability evaluated following 
the approach of lFullerton et al.l(ll996h . We clearly see the double 
peaked structure in the stellar lines. 



Previously estimated to 118 +25 km s"' (IConti & Ebbetsl 

11977 ) and 86kms ' (Penny, 1996), the projected rotational 
velocity of HD46150 has been determined using the Fourier 
method on four different lines (Hen/i4686, Oiii/15592, 



Civ /1 5801 and Hei/(5876). We find a mean value of v sin/ = 
97 ±9kms-i. 

In the literature, the spectral type is quoted b etween 
05 V ((f)) dUnderhill & Gilrov, 1990; Mafz-Apellani z et al 



2004 ) and 05.5V((f)) (Conti & Ebbets. ,1977; ,Garmanv et al. 



1980 ). From our measurements of the EWs, we confirm the latter 



spectral type. We observe weak Niii/l/l 4634-41 emission lines 
and a strong He ii A 4686 absorption line which confirms the ad- 
dition of the ((f)) suffix. We compute My - -5.45 + 0.22, from 
V - 6.74 and ^(B-y) - 0.40, which corresponds to a MS star 
although we note that My is s omewhat brighter than the theoret- 
ical value (M{J"=° = -5.12) of M artins & PlezI (l2006h for typical 
05.5 V stars. The difference could result from the approximative 
distance taken at 1.55 kpc or could also reflect the presence of a 
companion. We conclude that the spectral type of HD 46150 is 
05.5 V ((f)). Despite the lack of a clear periodicity or of a clear 
detection of the companion, several hints (double peaked struc- 
tures in TVS, large ARV,...) lead us to consider HD 46150 as a 
potential binary. 



3.3. Presumably single stars 
3.3.1. HD 46056 

In the past, HD 46056 was generally considered as a spec- 
troscopic binary because of a variable RV ( Walbornl 119731 : 
ILiu et all 119891: lUnderhiU & Gilrovi 11990 ). The fine profiles 
always appear broad and shallow, suggesting the star to be a 
single rapid rotator rather than a spectroscopic binary. Indeed, 
the Helium lines display very broad profiles extending over 
more than 15 A with a clearly non-Gaussian shape (see Fig.|2ll. 
Furthermore, some metallic lines, such as Mgii/14481 or 
OmA 5592, are too shallow to allow a measurement of their RVs 
by Gaussian fits or are too heavily blended to be distinguished as 
it is the case for the C iv AA 5801-12 lines. As a consequence, we 
adopted a different technique to measure the RVs. We first com- 
pute the Fourier transform of different Helium line profiles using 
our best spectrum of the star which combines a high signal-to- 
noise ratio and a high dispersion. We then found the mean pro- 
jected rotational velocity ( v sin /) to be 355 ± 2 1 km s" ' . 

Synthetic rotation profiles with v sin / = 346, 350, 358 and 
365 km s"' were subsequently used to fit the He 1/14471, 
He II A 4542, He ii A 4686 and He i A 5876 lines, respectively. For 
all these lines, the mean of the differences between the two ex- 
treme RV values were 28.1 kms"' and the mean 1-cr disper- 
sion is close to 9.5 km s"' . This is close to the adopted measure- 
ment error for such broad lines. However, the line profiles appear 
to vary significantly as confir med by the TVS a nalysis. This is 
reminiscent of HD 326331 (see lSana et a l.L 12008 ) and we there- 
fore conclude that the lines of HD 46056 display some variations 
though they are not easily associated with the presence of a com- 
panion. For instance, non-radial pulsations could be responsible 
for changes in the line profiles of HD 46056, likely to perturb the 
measurement of the radial velocities. More spectra and most of 
all a dedicated monitoring (spectroscopic as well as photomet- 
ric) with a high temporal resolution are necessary to probe the 
existence of such short-time variations attributable to non-radial 
pulsations. The detection of such variations in t he line profiles of 
rapid r otators has, for instance, been reported bv lDe Becker et all 

dloosh . 

To estimate the spectral classification, as the line profiles 
are clearly not Gaussian, we measured the EWs by integrat- 
ing the line profiles. Conti's criteria yield a spectral type 08 V. 
Moreover, the visual absolute magnitude of this object, calcu- 
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Fig. 7. He I /1 5876 line profiles of HD 46056 observed from 
February 2006 to March 2008. 



lated using V = 8.16 and E{B - V) = 0.47, is about My = 
-4.25 + 0.21 which corresponds to a MS star as suggested b y the 
EWs of the classification lines. Unlike lR'ost & Conti (1 197^ who 
noticed some sporadic emissions in the Ha line and assigned 
a spectral designation e for this star, we report that the narrow 
emissions, visible in the Balmer H/3 and Ha lines, are probably 
caused by the presence of nebular material around HD 46056. 
This conclusion is supported by the presence of other nebular 
lines, e.g. the [Om]/I5007 emission line. Therefore, a spec- 
tral type 08 V n i s assigned to HD 4 6056, in agreement with 
IWalbornI (UqtI or lUnderhiU & Gilrov> (.1990.) . 



3.3.2. HD 46202 

Alth ough some au t hors claim that HD 46202 may b e RV vari- 
able dAbt & Biggsl 119721: lUnderhill & Gilroyl Il990l) . we found 
no variation from our two year spectroscopic campaign (between 
January 2006 and January 2008) and we therefore conclude that 
HD 46202 is a presumably single star. Indeed, the RV disper- 
sions are similar to those measured for the interstellar lines and 
the average RVs agree with those reported in the literature. Our 
results therefore s upport the previous investi gatio ns made by 
iPlaskett & Feared (fl931h and lPetrie & Feared (Il96lb . 

The projected rotational velocity of HD 46202, estimated 
using the Fourier method, amounts on average to v sin/ - 
54 + 15 km s"' for the stellar lines. 

Previous studies reported HD 46202 as an 09 V star 
(Morgan etal .L I1955L Il965t IWalbornL I197U IConti&Leet^ 
[l974:,Massev et al.L 119951) ." Our data confirm this classification. 



Moreover, we obtained Mv = -4.13 + 0.21 from V = 8.18 and 
E(B - y) = 0.44, which indeed corresponds to a MS star. 

A s quoted in the literature (Per ez, 1991^, iMa ssev et al.L 

119951: iMaiz-Apellaniz et al.L |2004|) . the visual magnitude seems 
constant in the range 8.18-8.21, with the exception of 



lOgura & Ishidal (Il98ll) who derived a V magnitude of 8. 10. 



3.3.3. HD 46223 

HD 46223 is the hottest member of NGC 2244. The star was 
regularly obse rved for studies of the interstellar medium (e.g. 
IZagurvir2001h but the qu estion of the variabil it y of th is star arose 
already in the 1970's. 'Underbill & Gilrov' 099^) concluded 
from published data (Cruz-Gonzalez et al., 1974; Conti et all 
[19771: iGiesl Il987h and from their own observations that the star 
may vary. The hypoth esis of a s pectro scopic binary was sup- 
ported by the study of iLiu et aL I (Il989h who noticed a RV dis- 
persion of 10.3 km s ' for three observations. However, these RV 
changes co uld also be interpreted as motions in the atmosphere 
of the star (lUnderhill & Gilrovll 19901) . 

The dispersions of the RVs of the more intense lines such 
as He II /1 4542, Hen /1 4686 and He 1/15876 are similar to those 
of the DIB. These values are slightly smaller than those of the 
metallic lines and the He 1/14471 line. In any case, the disper- 
sion is never significant even for these metallic lines: indeed, 
the larger noise of these faint lines might explain their slightly 
higher RV dispersion. 

In consequence, we also applied the TVS method to our 
spectral time series. We performed this test with, on the one 
hand, only the Aurelie data, and, on the other hand, with all 
collected spectra (treated to get similar resolution, see Section 
13.2. It . In both cases, we failed to detect any significant vari- 
ations around the He 1/14471, /1 5876, Hen /1 4542 and near the 
C iV/l/15801-12 lines (Fig.|8]l and we therefore consider the star 
to be constant. 
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Fig. 8. Mean spectrum and TVS of HD 46223 computed from 
the Aurelie data for the He 1/14471 (on the left) and the 
He II /1 4542 lines (on the right). The dotted line illustrates the 
99% significance l evel for the vari ability evaluated following the 
approach of iFullerton et al.1 (119961) . 
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To determine the projected rotational velocity of the star, we 
used the Fourier transform on the same lines as for HD 46150. 
We obtain a mean projected rotational velocity of about v sin / - 
100 + 17 km s~' while the literature reports a value of v sin / - 
103kms-i (Pennvi 119961) . 



HD 46223 was previously classi fied as 05 ( Hiltneit 1956t 
iJohnsonll 1 962l:lBisiacchi et aU I 1 9821). 04 ([Morgan et all 1 19651) . 
05((r)) (IConti & Leed 119741) and, finally, as 04V((f)) by 
iMassev et al.1 ( Il995h . The measured EWs correspond to an 04 
star. Using V - 7.27 and E(B - V) - 0.50, we compute a vi- 
sual absolute magnitude of about My = -5.23 + 0.22, favoring 
a MS classification even though this va lue is slightly brighte r 
than the theoretical one (M{^™ = -5.56.l Martins & Plezll2006l) . 
Moreover, the spectrum clearly presents strong Nm/t/l 4634-41 
emission lines accompanied by strong He n A 4686 absorption. 
HD 46223 also displays both Si iv AA 4088 and 41 16 in emission. 
Therefore, we assign an 04 V ((f^)) spectral type. We also note a 
weak emission in the Ha line probably due to the nebular mate- 
rial surrounding the star as we already mentioned for HD 46056. 



3.3.4. HD 46485 

In previous works on NGC 2244, little attention was paid to 
HD 46485. The literature presents only a few measurements but 
in no way the results of a dedicated observing run. 

As for HD 46056, HD 46485 presents very broad and moder- 
ately deep lines suggesting that this star is a rapid rotator (Fig|9]|. 
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In order to determine the projected rotational velocity of 
HD 46485, we independently applied the Fourier transform on 
three different lines: Hen /1 4686, On: /1 5592, He U 5876 and 
we found a mean value of 301 +25kms"^ As the line profiles 
are not Gaussian, we used this measurement to generate sev- 
eral theoretical rotation profiles with v sin ; - 325, 3 1 8, 29 1 and 
301 km s ' which we correlate with the observed He 1/14471, 
He II /1 4542, Hen /1 4686 and He 1/15876 lines, respectively, to 
derive the RVs. The 1-cr range of these lines are similar to each 
other (see Table[3]l and, as for HD 46056, similar to the RV mea- 
surement error. However, obvious profile changes are clearly de- 
tectable though they are difficult to attribute to a possible binarity 
(Fig.|9l). They may be due to non-radial pulsations. 

IConti &Leep| (11974 identified HD 46485 as an 07.5 star 
while IWalbornI il9lW suggested an 08 V n(e) classification. 
Using Conti's criteria, we infer a spectral type of 08 V for 
HD 46485. With y = 8.20 and £(B - V) = 0.58, we derive 
a visual absolute magnitude M y - -4.68 + . 36. T his value 
agrees with a typical 08 V star dMartins & Plezl l2006h . In con- 
sequence, with its broad diffuse lines, we classify HD 46485 as 
an 08 V n star. 



3.3.5. HD 46966 

For our study about the variability of HD 46966, we collected 22 
spectra between January 2003 and January 2008. HD 46966 does 
not present any significant RV dispersion neither for the metal- 
lic nor for the Helium lines (/1 4471, /1 4542, /1 4686 and /1 5876). 
Moreover, no significant variation is observed by the TVS 
method. This suggests that HD 46966 is a presumably single star. 

iPennvl (11996*) estimated the projected rotational velocity of 
HD 46 966 close to 59kms '. Later on, Munari & Tomasell^ 
(Il999h revised the value to 90 km s ' . From four different lines 
(already mentioned for HD 46573), we measured the v sin ; with 
the Fourier method and the obtained mean value is estimated 
to be 63 + 17 kms~^ in good agreement with the value given by 
IPennvl (119961) . 

The spectral typ e s quo t ed in the literature are 8 V 
(iMunari & Tomasella'. '1999'; 'Mafz-Apellaniz et all |2004 or 
OB 5 V (Conti et al., 1977; Garmany et al., 1980!). Based on the 
Conti criteria, we derive an 08.5 V spectral type. Although lo- 
cated outside NGC 2244, we assume that HD 46966 is situated 
approximately at the same distance as the open cluster (1.4-1.7 
kpc). As before, we compute, with V - 6.87 and E(B - V) = 
0.21, the minimum visual absolute magnitude My = -4.51 and 
the maximum one at about My = -4.93. Thes e two values are 
larger than the theoretical one (M{J"=°= -4.27. iMartins & Pled 
l2006h . However, we consider this difference as an effect of the 
uncertain distance rather than a consequence of binarity and we 
propose a MS classification for HD 46966. 

3.3.6. HD 48279 
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Fig. 9. Same as FigQ but for HD 46485. 



Previo us spectroscopic investigations by iPetrie & Reared 
(1 19611) reported a r nean r adial velocity of 31.0 + 5.1kms"'. 
lUnderhill & Gilrovl (Il990l) . on the basis of their measured 
RV an d al so of those of iP etrie & Pearcd (Il96lh . IConti et al.] 
(1 19771) and lSanford & Merrill (.1938) , suspected the star to be a 
single-fined spectroscopic binary. 

We analyzed 4 spectra in the yellow domain and 10 in the 
blue region, along with 6 echelle spectra covering the whole 
visible domain, therefore including the blue and yellow wave- 
bands (see Table [U. The differences between RVs measured on 
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Helium lines are not large enough to support the binary scenario. 
Although the variations are much stronger in the fainter metallic 
lines than in the stronger He i A 447 1 and A 5876 lines, the accu- 
racy on the RV measurements of the more shallow lines is more 
affected by the noise. In addition, no signature of a companion, 
such as asymmetries in line profiles, has been found in any spec- 
trum of our time series. The TVS method, first applied to spectra 
in the blue domain and then to the data in the yellow region, did 
not show any significant variations for He 1/14471, He 1/15876 
nor for the metallic lines. The likelihood that HD 48279 is a bi- 
nary thus remains low. 

We estimated the projected rotational velocity from the 
Fourier method to v sin; = 130+ 13km s '. This value 
is in agreement wi th t hat of 123k ms~' determined by 
IConti& Ebbetsl(ll977h and lPennvl(ll996h . 

The spectrum of HD 48279 presents, beside the Balmer lines, 
deeper He i lines than He ii ones, the obvious signature of a late 
O-type star and we infer an 07.5 classification from Conti's cri- 
teria. We also find Si iV/1/1 4088-4116 in deep absorption. The 
most remarkable features of the spectrum are the prominent N iii 
absorption lines at 4195, 4634 and 4641 A, implying an ON 
spectral designation. 

IOgura&Ishidal(ll981h do not mention HD 48279 as a mem- 
ber o f NGC 2244 but it still belongs to the Mon OB2 associa- 
tion dde Wit et al.L .2004) . Without better information, we con- 
sider the star to be at the same distance as NGC 2244. Using 
V = 7.91 and E{B - V) = 0.407, we compute the two extreme 
visual absolute magnitudes for a distance between 1.4 and 1.7 
kpc. We obtain My = -4.08 and My - -4.50 respectively, sit- 
uating HD 48279 on the MS band. Accounting for the strength 
of the nitrogen lines, we hence adopt the ON7.5 V spectral type 
for HD 48279 slightly ear lier than previo us estimates (ON8 V, 
IConti & LeeD..1974: 08 V, IWalbornlll973"l) . 



3.4. Photometry and close "visual" companions 

In addition to our results, it is interesting to look at the photomet- 
ric variability of each star We thus retrieved, from the Hipparcos 
and the ASAS (All Sky Automated Survey) databases, the 
data corresponding to each studied object. These two datasets 
have been analyzed independently since the instrumentations are 
quite different. 

Table|5]lists the mean V magnitude of each star and the 1-cr 
standard deviation (on V) of the best data provided for each star 
in the ASAS catalog. The Fourier transform applied to each stud- 
ied object reveals no significant periodicity, concluding there- 
fore that photometric variations of the stars are absent on the 
observed time series. 



Table 5. The mean V magnitude and the l-cr standard deviation. 



Star 


Mean V magnitude 


l-cr error 


HD 46056 


8.132 


0.006 


HD 46149 


7.598 


0.006 


HD 46150 


6.720 


0.021 


HD 46202 


7.946 


0.048 


HD 46223 


7.270 


0.006 


HD 46485 


8.255 


0.011 


HD 46573 


7.925 


0.008 


HD 46966 


6.854 


0.019 


HD 48279 


7.596 


0.026 



The Hipparcos catalog provides us with a sufficient number 
of data for only four stars of our sample: HD 46150, HD 46485, 
HD 46573 and HD 46966. For these objects, we compute a mean 
magnitude and a corresponding l-cr standard deviation on the 
magnitude of 6.787 + 0.007, 8.342 + 0.023, 8.030 + 0.011 and 
6.851 + 0.007, respectively. Once again, these small dispersions 
do not allow us to conclude to the existence of any significant 
variability for any of studied stars. 

The speckle interferometric campaign of Galactic O-type 
stars with V < 8 undertaken by iMason et al.l (Il998 h revealed 
that HD 46150 and HD 48279 have "visual" companions whose 
separations range from 2.7"-74.6" and 4.6"-56.2", respec- 
tively. However, we note that the "visual" companioi i close to 
HD 48 279 was quoted as optical and not physical by iLindroosI 
(^19 851). In ad dition, the Washington Visual Double Star (WDS, 
Mas on et"al1 12001) catalog also reported a "visual" compan- 
ion located between 9.6" and 10.5" of HD 46056. However we 
stress that no evidence exists so far to confirm with any cer- 
tainty the physical association of these companions, especially 
with such high angular sepa rations. 

More recently. Turner et al.l ( l2008l) published results of an 
adaptative optic survey to search for faint companions among 
Galactic O-star systems (with V < 8). They observed in the 1- 
Band in order to detect companions in the projected separation 
between 0.5" and 5.0". They observed five stars of our sample 
(HD46149, HD46150, HD46223, HD46966 and HD48279) 
but found no companion. 

4. Discussion 

4.1. Multiplicity of O-type stars 

In this paper, we have studied all six O-stars of NGC 2244: 
HD46056, HD46149, HD46150, HD46202, HD46223 and 
HD 46485. HD 46149 is clearly a long-period SB2 while 
HD 46150 is a good binary candidate. The minimum frequency 
of spectroscopic binaries among O-stars in this cluster is thus at 
least 17% and could possibly reach 33%. 

4.1 .1 . The detection biases 

Because of the spectroscopic approach, a bias exists towards the 
detection of O+OB binary systems, providing only a lower limit 
on the true binary fraction. Indeed, several observational biases 
could make us miss a binary system with either a low orbital 
inclination, a rather large mass and/or brightness ratio or a very 
wide separation (and though either very long orbital period or a 
large eccentricity). 

We first use the approach of i Garmanv et alj ( |1980|) to esti- 
mate the probability that we have missed a binary system be- 
cause of a low orbital inclination as a function of the assumed 
mass ratio and orbital period. In this scheme, we consider that 
the semi-amplitude of the RV curve K should be smaller or equal 
to twice the RV dispersion (cr^y) of the photospheric lines, i.e. 
He II /1 4542 in the blue domain and C iv/l5812 in the yellow re- 
gion (or He I /1 5876 when this metallic line is not visible). The 
mass function of the primary component is given by: 



f(m) = 



Ml sin i 

1.0355 X 10"^ P(l 



(1) 

(2) 



where Mi is the primary mass (in Mq), K is expressed in km s ' , 
e represents the orbital eccentricity, q - M1/M2 and P defines 
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the orbital period (in days). 

Equation |2] is transformed to express the orbital inclination as a 
function of other parameters. To obtain an upper limit on sin /, 
we insert Ictrv as an upper limit on K and assume a zero eccen- 
tricity. Therefore, we obtain: 

sin/ < 9.392 x 10^^ tr«v (^^"^g^)'^' (3) 

Making the assumption of a random distribution of the orbital 
directions in space, we can write the probability that the orbital 
inclination is smaller than the value obtained from Equation [3] 
as: 

nil,,, 

I sin it// = l-cos;„p. (4) 
Jo 

The probabilities reported in Table 00 are the mean values com- 
puted from the two photospheric line s quoted above and by using 
as primary mass, the masses listed in lMartins et alj (l2005h . 

The probabiHties that HD 46202, HD 46966 and HD 48279 
{^2QMq) would be binaries with an orbital period less than 
28 days and a companion earlier than A5 {-2Mq, q- 10) are 
very low (i.e. < 7%). Any missed binary would most probably 
have to be a system with a rather large mass ratio and a long or- 
bital period. 

The same conclusion could be reached for HD 46223 (^47 Mq). 
At a period of 28 days, the probability of having an O-type com- 
panion is close to zero (~ 0.2%) and we have only 3% chance 
to miss a B3 secondary (^ 8 Mq, q - 5). However, for the same 
orbital period, there is a 10% chance that the binary signature 
could have been missed if the spectral type of the second com- 
ponent would be B6 or later (< 5 Mq, q- 10). 

We also perfomed Monte-Carlo simulations following the 
method of Sana et al. (2009, submitted). The orbital parameters 
are randomly generated ( 1 0000 trials) under the assumptions that 
50% of the binaries have an orbital period of P < 10 days while 
the other 50% have f > 10 days. Moreover, the periods are uni- 
formly drawn in log P space (between 0.3 and 3.5), the system 
orientations (cos /) are randomly drawn among a uniform distri- 
bution, the mass ratio q-M\ /M2 is also uniformly distributed 
between 0.1 and 1.0 and we uniformly pick up eccentricities 
between 0.0 and 0.8. The simulations show that the probabil- 
ity (Table 1^ to have missed, in HD 46202 and HD 46223, a 
peak-to-peak RV variation (ARV) larger than the observed value 
of 20km s"' is less than 1% for systems with 2<P< 10 days 
and less than 10% for 10 < f < 365 days. For the rapid rotators, 
the simulations were run for ARV > 40 km s ' since the measure- 
ment errors are larger The probabilities amount to 2% and 20% 
for HD 46056 and HD 46485, respectively. Therefore, it is very 
unlikely that short-period O-i-OB systems were missed by our 
study, and this result does not depend on the small number of 
stars in NGC 2244. For long-period (365-3000 days) systems, 
the probability of missing binaries increases to 50% (70% for 
rapid rotators). This is expected since our sampling does not per- 
mit us to investigate thoroughly this part of the parameter space. 

' Table |6] reports the probabilities that a binary system could have 
been missed for different values of the orbital period P and of the mass 
ratio q= M1/M2. 

Table|7]lists theprobability of missing a significant orbital RV vari- 
ation (ARV > 40 km s"' for the two rapid rotators and ARV > 20 km s"' 
for the other stars). Mi represents the primary mass and is expressed in 
Mq. The last four columns are classified according to the orbital period 
of the system (in days). 



Table 7. Probabihty of missing a significant orbital RV variation. 



Star 


Ml 


2-10 


10-365 


365-3000 


2-3000 


HD 46056 


20.8 


0.02 


0.17 


0.77 


0.24 


tiU 46149 


on 
20.0 


O.Oi 


0.04 


U.3I 


U.U/ 


HD 46150 


34.4 


0.01 


0.03 


0.21 


0.05 


HD 46202 


17.1 


0.01 


0.09 


0.50 


0.13 


HD 46223 


46.9 


0.01 


0.03 


0.22 


0.05 


HD 46485 


20.8 


0.01 


0.20 


0.68 


0.20 


HD 46573 


22.9 


0.01 


0.06 


0.38 


0.09 


HD 46966 


18.8 


0.01 


0.06 


0.34 


0.08 


HD 48279 


22.9 


0.01 


0.07 


0.40 


0.10 



4.1 .2. The effects of stellar evolution 

We now compare our results with those of other young open 
clusters studied in a similar way. The minimal binary fraction 
in NGC 2244 (17%) appears close to that of IC 1805 (20%). 
Howeve r, we note differences in the stellar content of these two 
clusters. [De Becker et al.l (1200 6*) have shown that at least one gi- 
ant (0 111(f)) and one supergiant (Olf) belong to IC 1805 while 
all O-type stars in NGC 2244 are MS objects. Evolution effects 
(and thus the age of the systems) therefore seem to play no role 
in the binary properties though to be firmly established, this con- 
clusion must be confirmed by studies of other clusters. 

Garcia & Mermilliod (2001) asserted that the orbital periods 
of spectroscopic binaries in the O-type star rich clusters were 
concentrated around 4-5 days. We emphasize that the absence 
of short-period binaries in NGC 2244 does not support this con- 
clusion (see previous subsection). 

A comparison between NGC 2244, IC 1805 and NGC 6231 (all 
three being now well-studied) is provided in Table[8] The first 
and second columns give the name and the density of the cluster 
The third column provides the minimal binary fraction deter- 
mined after intensive monitoring of the O-type stars. The next 
columns yield the number of 0-stars in each cluster, the number 
of detected short and long-period binaries and the variable stars 
(potential binaries and probable intrinsic variables). The last col- 
umn quotes relevant references. Differences are readily seen in 
this table. Apart from the already quoted difference in the pro- 
portion of short vs. long-period binaries, there is also a differ- 
ence in the preferential nature of the companion. If we consider 
all O-type stars to be binaries, those apparently single having 
undetectable low-mass companions, the majority of O-stars in 
NGC 623 1 appears to have OB companions, while the signature 
of additional early-type objects is found in only a few O-stars of 
NGC 2244 and IC1805. 

4.1 .3. A possible impact of the stellar density 

The surface density of a given cluster is evaluated through the 
ratio between the number of members of the cluster and the cir- 
cular area del ineated by a typical radius, e.g. the one given by 
iTadross et ai](i2002,) . Both IC 1805 and NGC 2244 have a rather 
loose cluster core whilst the stellar surface density appears much 
higher in the central part of NGC 6231. With the new binary 
fraction from De Becker et al. (2006), Sana et al. (2008) and this 
work, the anti-correlation between the c luster density and its 
massive star binary fraction proposed by Garcia & MermilliodI 
(1200 Ih can no longer be supported. On the contrary, a possible 
correlation might better reflect the behaviour detected in these 
three clusters. However, we must remain careful since this as- 
sumption is directly linked to the error bars on the binary frac- 
tion. If it is confirmed, this would emphasize the important role 
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Table 6. Probability to have missed a binary system. 



q=\ 



P = 5 days 

^=5 q=W 



P=\A days 
q=\ q=5 q=lO 



P = 2S days 
= 1 ^=5 ^=10 



HD 46056 3.9E-3 5.0E-2 19.4E-2 

HD 46150 3.6E-3 4.6E-2 18.1E-2 

HD 46202 4.7E-4 6.0E-3 2.2E-2 

HD 46223 6.4E-4 8.1E-3 2.9E-2 

HD46485 5.1E-3 6.7E-2 26.6E-2 

HD 46966 3.4E-4 4.3E-3 1.6E-2 

HD 48279 2.3E-4 2.9E-3 l.OE-2 



7.8E-3 10.3E-2 44.9E-2 

7.1E-3 9.5E-2 45.5E-2 

9.4E-4 1.2E-2 4.3E-2 

1.3E-3 1.6E-2 6.0E-2 

l.OE-2 13.8E-2 71.5E-1 

6.8E-4 8.6E-3 3.1E-2 

4.5E-4 5.7E-3 2.0E-2 



1.2E-2 16.9E-2 

l.lE-2 15.7E-2 34.7E-2 

1.5E-3 1.9E-2 7.0E-2 

2.0E-3 2.6E-2 9.7E-2 

1.6E-2 23.1E-2 

l.lE-3 1.4E-2 5.0E-2 

7.1E-4 9.0E-3 3.3E-2 







^ = 60 days 






P = 300 days 






q=l 


q=5 


^=10 


q=l 


q=5 


q=lO 


HD 46056 


2.1E-2 


30.4E-2 




6.1E-2 






HD 46150 


1.9E-2 


28.8E-2 


78.2E-2 


5.7E-2 


53.3E-2 




HD 46202 


2.5E-3 


3.2E-2 


12.1E-2 


7.3E-3 


9.8E-2 


51.5E-2 


HD 46223 


3.4E-3 


4.4E-2 


17.2E-2 


9.9E-3 


13.7E-2 


19.6E-2 


HD 46485 


2.7E-2 


43.3E-2 




8.2E-3 






HD 46966 


1.8E-3 


2.3E-2 


8.5E-2 


5.3E-3 


6.9E-2 


29.0E-2 


HD 48279 


1.2E-3 


1.5E-2 


5.5E-2 


3.5E-3 


4.5E-2 


17.1E-2 



Table 8. Summary of the multiplicity in young open clusters. 



Cluster 


Density 
(stars/pc^) 


Minimal binary 
fraction 


n# O-stars 


Short-period 
binaries (P < 10) 


Long-period 
binaries 


Variable 
stars 


References 


NGC 2244 


2.84 


17% 


6 





1 


3 


This paper 


IC 1805 


1.65 


20% 


10 


1 


1 


4* 


Rauw & De Becker (2004) 
De Becker et aL (2006) 


NGC 6231 


12.46 


63% 


16 


6** 


4 


3 


Sana et al. (2008) 



* Notes: Four late 0-type stars were not studied by the two papers mentioned in the references. We will consider these objects as variable in the 

absence of a thorough study. 
'* Notes: We take into account that NGC 6231 also contains a WR+O binary since a WR is an evolved O-star. 



of dynamical gravitational inte ractions in the formation of mas- 
sive binaries ( IZinneckeill2003l) . To be ascertained, however, this 
relation needs to be confirmed by more data. The extreme case 
of loose OB associations might provide a good testbed to this 
aim. 

4.2. Rotational velocities 

HD 46056 and HD 46485 are two rapid rotators whose lines ap- 
pear slightly variable. Rapid rotators could either 

• be born as such and maintain a rapid rotation over their main- 
sequence lifetime, 

• be spun-up through transfer of mass and angular momentum 
during a Roche lobe overflow in a binary system, 

• be spun-up due to the core contraction at the end of the main- 
sequence. 

For these two rapid rotators, it seems more likely that the first 
hypothesis is the correct one since we found no clear indication 
for multiplicity and there is no evidence that these objects are 
at or near the end of their main-sequence lifetime (chemical en- 
richment,...). 

We also compared our distribution of the projected rotational 
velocities of the O-type stars in and around NGC 2244 to the cor- 
responding distribution for B-type stars. Indeed. iHuang & GiesI 
(2006) performed a survey of projected rotational velocities for 
a large sample of mainly B-type (later than 09.5) main-sequence 
or giant stars in young open clusters, including NGC 2244. The 



V sin / values were derived through comparison between theoret- 
ical line profiles and observe d profiles of the He i AA 4026, 4387, 
4471 and Mgn/(4481 lines. iHuang & GlesI (l2006h report pro- 
jected rotational velocities of 41 stars in NGC 2244. The aver- 
age projected rotational velocity of the B-type stars amounts to 
168 kms '. The 15 stars that are earlier than B3 show an aver- 
age V sin/ of 153 kms"'. Five stars out of the total sample (i.e. 
12 %) were found to display a projected rotational velocity above 
300 km s"' . At first sight, the distribution of rotational velocities 
of the O-type stars looks diff'erent from that of the B-type stars 
(see Fig. [Toll. We have performed a two-sample Kolmogorov- 
Smirnov test to check whether or not the two distributions are 
significantly different. On the basis of this test and adopting a 
usual significance level of 0.05, the null hypothesis (both sam- 
ples are drawn from the same parent distribution) cannot be re- 
jected. This is true for both the full sample and the six O-stars of 
NGC 2244. 

5. Summary and conclusions 

Our long-term monitoring has shed new light on the multiplicity 
of six massive O-stars in NGC 2244 plus three in the association 
Mon OB2. Table ISlFI summarizes the present knowledge on these 
stars. 

^ The first column provides the name of each star studied in this pa- 
per. The next two columns give the spectral classification and the pro- 
jected rotational velocity (in kms"') of the O-type stars while the last 
column reports the status of the star as derived from our campaign: "C" 
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0.4 

0.3 

0.2 

0.1 


0.4 

I 0.3 

^ 0.2 

(U 

£ 0.1 


0.4 



0.3 
0.2 
0.1 - 






41 B stars 
(Huang & Gies) 



2b single B stars 
(Huang & Gies) 



9 — stars 
(this paper) 



Table 9. Summary of the optical properties of O-stars 
NGC 2244 (top part) and Mon OB2 (bottom). 



in 



100 200 

V sin(i) 




400 



Fig, 10. Top : the distribution of projected rotational velocities of 
(1) the full sample of B stars studied by Huang & Gies (2006, top 
panel), (2) those B stars of Huang & Gies that show a constant 
radial velocity (middle panel) and (3) our sample of O-type stars. 
Bottom : the cumulative distributions of the three samples: all B- 
stars (solid line), single B stars (dashed-dotted line) and O-stars 
(dashed line). 



In the association Mon OB2, the cuiTent study reveals the 
detection of a new SBl: HD 46573. Although we determine a 
period of 10.67 days and we compute the orbital parameters of 
the system, these results must still be considered as preliminary. 
Indeed, the accuracy of the RV measurements could strongly in- 
fluence the orbital solution of the system, considering the low 
amplitude of the radial velocity variations. The question of the 
eccentricity of the orbit remains open but we note that the eccen- 
tric radial velocity curve appears slightly better than the circular 
one. 

Another interesting result, in NGC 2244, is the detection, for 
the first time, of a secondary component in HD 46149. However, 
for the time being, we can only constrain the period to be at 



means that our study did not reveal any indication of RV variability 
or of binarity, suggesting the star is probably single, "Bin?" indicates 
a potential binary, "SBl" corresponds to a spectroscopic binary where 
the second component has not been detected yet but where a period has 
been determined, whilst "SB2" corresponds to a spectroscopic binary 
where the two components are visible. 



Name 


Spectral Type 


V sin/ 


Spectroscopic status 


HD 46056 


OSVn 


355+21 


C ? (rapid rotator) 


HD46149 


08V+BO-1V 


78± IP 


SB2 


HD46150 


05.5V((f)) 


97 ±9 


Bin? 


HD 46202 


09V 


54 ±15 


C 


HD 46223 


04V((f+)) 


100 ± 17 


C 


HD 46485 


OSVn 


301+25 


C ? (rapid rotator) 


HD 46573 


07.5V((f)) 


110 + 18 


SBl 


HD 46966 


08.5V 


63 + 17 


C 


HD 48279 


ON7.5V 


130 + 13 


C 



* Notes: The reported value for HD 46149 is the projected rotational 
velocity for the primary. 



least several tens of days. HD 46149 is considered, by virtue of 
its mass and flux ratios, to be an 0+B binary. An observation 
campaign should be initiated to monitor the star and better con- 
strain the orbital solution of the system. Finally, we also detect 
weak variations in the RVs for one star belonging to NGC 2244 
(HD46150) and variations of the line profiles, most probably 
unrelated to binarity, in the two rapid rotators (HD 46056 and 
HD 46485). The observed binary fraction is thus 17% at mini- 
mum. 

This work, combined to the results of studies of NGC 623 1 
and IC 1805, also revealed significant differences between clus- 
ters regarding the preference of O+OB systems and the distribu- 
tion of short vs. long-period binaries. These differences can not 
be explained by the cluster age but could be linked to density 
effects. Of course, to confirm such a general conclusion, it is im- 
portant to rely upon a large sample of well-studied very young 
open clusters that harbour a significant population of early-type 
stars. To date however, only a few clusters have been subject 
to such an intensive search for binary systems. It is therefore 
of paramount importance to perform similar studies on a larger 
sample of clusters, spanning as wide a range of physical proper- 
ties as possible. 
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Table 1. Observations used for the study of massive stars in both NGC 2244 and Mon 0B2. 



Object 


Cluster 


Obs. ran 


Telescope/Instrament 


Wavelength domain 


n# Spectra 


AT 


HD 46056 


NGC 2244 


Feb. 2006 


OHP- AUREL1E/I.52m 




5500-5900A 


6 


5.02 






Apr. 2007 


SPM - ESPRESSO/2.20m 


4450-*900A 


5 


6.00 






Nov. 2007 


OHP - AURELIE/1.52m 


1 








Jan. 2008 


OHP- AUREL1E/I.52m 


5500-5900A 


4 


3.00 






Mar. 2008 


SPM - ESPRESSO/2.20m 




1 




HD 46149 


NGC 2244 


Dec. 1994 


Asiago - Echelle/1.82m 




1 








Dec. 1998 


Asiago - Echelle/1.82m 




1 








Jan. 2003 


OHP - EL0D1E/I.93m 




1 








Nov. 2005 


OHP - ELODlE/1.93m 




1 








Jan. zuuo 


JioU - rHKUo/Z.ZUm 




1 








reD. ZUUO 


OXJT> ATTDT3T TT3/1 <0»^ 


jjuu— jyuuA 





< 1 A 
J.IU 






Anr ?007 


SPM - FSPRFSSO/7 70m 


CCAA CAAA A 

5jUU-j9UUA 


'X 

J 


3 97 






Apr. ZOO / 


OHP - AUKbElE/l.Dzm 


1 

1 








Nov. 2007 


OHP - AURELIE/1.52m 


4450^900A 


6 


25.05 






Jan. 2008 


OHP- AUREL1E/I.52m 


5500-5900A 


4 


5.03 






Mar. 2008 


SPM - ESPRESSO/2.20m 




5 


4.93 






Sep. 2008 


OHP - AURELIE/1.52m 


4450^900A 


1 








Oct. 2008 


OHP- AUREL1E/I.52m 


4450-*900A 


3 


4.02 


HD 46150 


NGC 2244 


Nov. 1999 


OHP-ELODlE/1.93m 




1 








Jan 2003 


OHP-ELODIE/1.93m 




1 








Nov 2004 


OHP - ELODIE/1.93m 




1 








Jan. 2006 


ESO - FEROS/2.20m 




1 








reu. ZUUO 


rl/oU - riiKUo/Z.ZUm 




1 








Feb. 2006 


OHP - AURELIE/1.52m 


5500-5900A 


6 


5.10 






Apr. ZUU/ 


orJVl - JiorKJiooU/Z.ZUm 




c 

J 


< QQ 






INOV. ZUU/ 


OTJD ATTDT^T TT^/I 


AA<f\ /IQAAA 
44jU— 4yUUA 


c 

J 


LJX} 1 






Tori OAAQ 

Jan. zuuo 


OTJT> ATTT>TirT TTI?/1 ^Ot^ 


^^AA ^QAAA 

jjuu— jyuuA 


4 


O QQ 

z.yy 






Mar 2008 


SPM - ESPRESSO/2 20m 




5 


3.95 






/-\„f OAAQ 

Uct. ZUUa 


rMJV ATTDT7T TT7/1 

Urlr - AUKriLlii/i.jZm 


AAZf\ ylOAAA 

44DU— 4yUUA 


A 

4 


AA 

4.UU 


riL) ^bZOZ 


JNLrC ZZ44 


Jan. ZUUo 


lioU - rtKUo/Z.ZUm 




1 
1 








Feb. 2006 


OHP - AURELIE/1.52m 


5500-5900A 


2 


1.98 






Apr 2007 


OHP - AURELIE/1.52m 


5500-5900A 


2 


5.00 






Nov. 2007 


OHP - AURELIE/1.52m 


4450^900A 


4 


20.00 






Jan. 2008 


OHP - AURELIE/1.52m 


5500-5900A 


4 


3.01 


HD 46223 


NGC 2244 


Jan. 1999 


OHP - ELODIE/1.93m 




1 








Nov 2005 


OHP - ELODlE/1.93m 




1 








Jan. ZUUo 


TI7C/~\ TI?TI?T>/~\C /O OAm 

iloU - riiKUo/Z.ZUm 




1 








Feb ?006 


FSO - FFROS/7 70m 












reb. ZUUo 


UHF - AUKbElli/l.jzm 


CCAA CAAA A 

jjUU-jVUUA 


6 


C A-^ 

J.U3 






Apr. 2007 


OHP - AlJRELIE/1.52m 


5500-5900A 


3 


6.00 






Nov. 2007 


OHP - AURELIE/1.52m 


4450^900A 


5 


21.98 






Jan. 2008 


OHP - AUREL1E/I.52m 


5500-5900A 


4 


3.06 






■\ /r OAAO 

Mar. 2008 


SPM - ESPRESSO/z.zOm 




3 


O Al 

2.01 






Oct. 2008 


OHP- AURELIE/1.52m 


4450^900A 


1 




HD 46485 


NGC 2244 


Nov. 2005 


OHP-ELODlE/1.93m 




1 








Jan. 2006 


ESO - FEROS/2.20m 




1 








Nov. 2007 


OHP - AURELIE/1.52m 


4450^900A 


6 


14.99 






Jan. 2008 


OHP - AURELIE/1.52m 


5500-5900A 


4 


3.05 






Mar. 2008 


SPM - ESPRESSO/z.zOm 




5 


4.02 


Tjpv /I 

liU 403 / J 


ivion yjDZ 


Ton OMMA 

Jan. ZUUO 


liovj - rrl/KL/o/Z.ZUm 




1 








reu. ZUUO 


iloU - rrLivUo/Z.ZUni 




1 








Feb. 2006 


OHP- AURELIE/1.52m 


5500-5900A 


3 


2.97 






Apr. 2007 


SPM - ESPRESSO/2.20m 


4450-4900A 


5 


5.99 






Nov. 2007 


OHP - AURELIE/1.52m 


6 


25.98 






Jan. 2008 


OHP- AURELIE/1.52m 


5500-5900A 


5 


3.05 






Oct. 2008 


OHP- AURELIE/1.52m 


4450-4900A 


3 


3.00 


HD 46966 


Mon OB2 


Jan. 2003 


OHP-ELODIE/1.93m 




2 


1.76 
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Nov 2004 


OHP-ELODIE/1.93m 




1 




Jan. 2006 


ESO - FEROS/2.20m 




2 


1.94 


Feb. 2006 


OHP - AURELIE/1.52m 


5500-5900A 


6 


5.01 


Feb. 2007 


OHP- AURELIE/1.52m 


5500-5900A 


1 




Apr. 2007 


OHP- AURELIE/1.52m 


5500-5900A 


1 




Nov. 2007 


OHP- AURELIE/1.52m 


4450^900A 


5 


24.05 


Jan. 2008 


OHP- AURELIE/1.52m 


5500-5900A 


4 


3.02 


UD 48279 MonOB2 Jan. 2000 


OHP-ELODIE/1.93m 




1 




Nov 2004 


OHP - ELODIE/1.93m 




1 




Nov. 2007 


OHP- AURELIE/1.52m 


4450-4900A 


9 


16.02 


Jan. 2008 


OHP- AURELIE/1.52m 


5500-5900A 


4 


3.01 


Mar. 2008 


SPM - ESPRESSO/2.20m 




4 


4.01 


Oct. 2008 


OHP - AURELIE/1.52m 


4450^900A 


1 





